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Abstract. - Recent experimental and theoretical studies have highlighted the possible role of a 
electronic nematic liquid in underdoped cuprate superconductors. We calculate, within a model 
of d-wave superconductor with Hubbard correlations, the spin susceptibility in the case of a small 
explicitly broken rotational symmetry of the underlying lattice. We then exhibit how the induced 
spin response asymmetry is strongly enhanced by correlations as one approaches the instability 
to stripe order. In the disorder-induced stripe phase, impurities become spin nematogens with 
a C2 symmetric impurity resonance state, and the disorder-averaged spin susceptibility remains 
only C2 symmetric at low energies, similar to recent data from neutron scattering experiments on 
underdoped YBCO. 



Introduction. — Incommensurate one dimensional 
composite spin and charge density waves, often called 
"stripes" PUS], have been observed and play an impor- 
tant role in discussions of the underdoped cuprates and 
other systems. Such states were predicted theoretically 
in the context of mean-field studies of Hubbard mod- 
els [TUTU], and later observed in neutron scattering ex- 
periments in La2- x Baa;Cu04 and L^-x-yNdySr-rCuC^ 
[TTinP2] . Stripes break the discrete translation and rotation 
symmetries of the Cu02 planes. Rotational symmetry is 
also broken in so-called liquid crystal analogs called "elec- 
tronic nematic" states, but these preserve translational 
symmetry and may occur as the initial instability of a 
paramagnetic state before an ordered state of charge, spin, 
or combined spin and charge order is reached [3| fT3l - fT6] . 

While stripe-like ground states were originally thought 
to be a very special feature of the 214 compounds, this 
view has changed in recent years, in particular with the 
discovery of broken C4 symmetry in the spin response of 
highly underdoped samples of YBa2Cu306+5 [E], anc b 
more recently, a subtle charge order in the same sys- 
tem [18| . Signatures of nematic order have also been 
recently reported in transport and tunneling measure- 
ments [T9H21] . In all these cases, however, the ques- 



tion of C4 — > C2 symmetry breaking is muddied by the 
fact that the crystal is not tetragonal, since for exam- 
ple in YBa2Cu3 06+5 the CuO chains give a well-defined 
anisotropy in the untwinned samples on which the ex- 
periments were performed, such that the system is for- 
mally orthorhombic. Nevertheless, the evolution from op- 
timally to highly underdoped samples, which is accom- 
panied by a dramatic enhancement of the anisotropy of 
the responses, is quite striking, and leads to the com- 
mon assumption that these highly correlated underdoped 
materials display a strongly enhanced "nematic suscep- 
tibility" , i.e. a tendency to create nematic order which 
is driven by the very small symmetry-breaking field pro- 
vided by the x — y anisotropy in the band-structure. How- 
ever, these ideas have rarely been cast in a concrete mi- 
croscopic model allowing a direct study of how disorder 
and local electronic correlations which drive the Mott in- 
sulating state influence nematicity. Previous studies have 
mainly focussed on the high-energy spin fluctuations of 
the RPA susceptibility of a homogeneous d-wave super- 
conductor with an anisotropic band-structure [22H26] . or 
utilized phenomenological Ginzburg-Landau approaches 
[iniHTlUB] . More recently, the nematic response has also 
been studied within the two-dimensional Hubbard model 
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with slight x — y hopping asymmetry using strong-coupling 
cluster methods, and found to be significantly enhanced by 
interactions at low temperatures in the underdoped pseu- 
dogap regime [2"W3"0] . 

The salient features of the inelastic neutron scattering 
experiments on strongly underdoped untwinned YBCO 
samples which should be reproduced by a reasonably com- 
plete theoretical analysis are as follows: 

1. The neutron intensity near (it, it) evolves from a pat- 
tern of four symmetrically placed incommensurate 
peaks at high energies which merge at (it, tt) at the 
spin resonance energy Qq, to a pattern with Ci sym- 
metry with two incommensurate peaks at quite low 
energies uj <C ^o- 

2. The details of this pattern appear to be important: 
the intensity in the nematic case has a saddle point 
form in q space, with a maximum at (tt, tt) for cuts 
along b and a minimum there for cuts along a. 

3. The strength of the C4 symmetry breaking increases 
as one underdopes. 

Some phenomenologies have been rather successful in ac- 
counting for some of these features, but to our knowl- 
edge there is no microscopic approach which has thus far 
successfully accounted for all of them. In this work we 
present a simple version of such a theory, based on our 
earlier work exploring the effect of Hubbard correlations 
on the d-wave superconducting state, to which we now 
add a small symmetry-breaking field to simulate e.g. the 
effect of the chains in YBa2Cu306+a. We show that the 
nematicity in the spin response is enhanced by correlations 
and decreasing temperature; what is perhaps more surpris- 
ing is our finding that the phenomenon can be enhanced 
further by pairing and disorder. We find that disorder is 
crucial to explain the saddle point structure of the inelas- 
tic scattering intensity in g-space. In addition, we present 
local investigations which exhibit the explicit formation 
of nematogens, nematically driven impurity states, which 
have not to our knowledge been observed in the cuprates, 
although recently reported in the stripe-ordered phase of 
Fe-based systems [3"TH3"3"] . Indications of an enhanced ne- 
matic susceptibility above the magnetic transition have 
also been observed in these systems [34] . 

Model. — The Hamiltonian is given by 
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+ Yl ( A *4^l+5l + AJiCi+i^Cit) , (1) 



where c ia creates an electron on site i with spin a, and 
tij = {t x ,t y ,t'} denote the hopping integrals to the two 
nearest neighbors. In Eq.([T]), and mi refer to the charge 
density and magnetization, respectively, Vi is an impurity 



potential from a set of N point-like scatterers, fj, is the 
chemical potential and is the d-w&ve pairing potential 
between sites i and j. The amplitude of is set by 
the superconducting coupling constant g and will exhibit 
a slight anisotropy inherited from a finite <5o = (t y —t x ). 
Below we fix the parameters t' = — 0.35i, adjust (i to give 
a hole doping x = 1 — 11 ~ 10%, and g = 0.6 leading 
to realistic pairing amplitudes A = O.lOi. The hopping 
asymmetry Sq = 0.05 is fixed for all results presented here, 
and we use units where t = t y = 1.0. We have solved 
Eq. (TT]) self-consistently on unrestricted N x N lattices by 
diagonalizing the associated Bogoliubov-de Gennes (BdG) 
equations at T = O.Oli. [35] 

The model given by Eq.([T]) has been used extensively 
in the literature to study the competition between bulk 
superconducting and magnetic phases, and field-induced 
magnetization [35] . It has also been used to study moment 
formation around nonmagnetic impurities in correlated d- 
wave superconductors [37] • In the case of many impu- 
rities, Eq.([T]) was used to model static disorder- induced 
magnetic droplets [55H^2] , and explain how these may 
increase in volume fraction when moving to lower dop- 
ing levels and eventually form a quasi-long-range ordered 
magnetic stripe phase. More recently, Eq.([T]) extended to 
the vortex state was used to obtain a semi-quantitative 
description of the temperature dependence of the elastic 
neutron response in underdoped LSCO [431144] . 

The transverse bare spin susceptibility Xo^O^i ^jj £a = 
—i J °° dte luJt ([erf (t),erj(0)]), can be expressed in terms 
of the BdG eigenvalues E n and eigenvectors u n ,v n as 
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Including the electronic interactions within RPA we find 
for the full susceptibility 

f»=2 [1 - UxT(«)]?t ?l xT{rur 5 ,uj). (5) 

Fourier transforming with respect to the relative coordi- 
nate r = Ti — fj defines the spatially resolved momentum- 
dependent susceptibility \ (q. R, w) = 53 -e"'' r ,Y(i?. f. w). 
Averaging over the center of mass coordinate R = (ft + 
fj)/2, this expression gives the susceptibility %(q, uj) rele- 
vant for comparison with neutron measurements. 

Results. — Figure Q] shows the susceptibility 
x"(q,uj) = Imx(q, oj) at low energy uj/t = 0.0 versus q x 
and q y for different U in the homogeneous superconduct- 
ing phase. As seen, the apparent asymmetry in the neu- 
tron response is clearly enhanced as the correlations are 
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Fig. 1: (Color online) Constant energy cuts of the susceptibility 
X"(g,w) at low energy u/t = 0.0 for S = 0.05, N = 80, and 
U/t = 2.0 (a), U/t = 2.2 (b), U/t = 2.4 (c), and U/t = 2.6 (d). 

0.4 ■ ■ ■ 1 
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Fig. 2: Nematic spin response r\ versus U for <5o = 0.05 in the 
law-U homogeneous d-wave superconducting phase. 

enhanced. Quantitatively, we can define a nematic spin- 
response "order parameter" as 

= max{x"(gx, 7T, u; )} - max{x"(7T, u; )} 
max{x"(g a ;, 7r, uj )} + max{x"(7r, q y , uj )} ' 

Figure [2] shows the dramatic increase of r] as U approaches 
the stripe instability from below for fixed hopping asym- 
metry 8q = 0.05. In cuprates with weak orthorhombicity, 
it is therefore natural to expect a significant enhancement 
of the nematic response upon lowering the doping level 
from the optimally doped regime fl"6l[29l l30 ] . 

Experimentally, it is known that the observed x — y 
asymmetry disappears as one increases the energy ui which 
seems natural in light of the presumably very small energy 
scale associated with the nematic aligning field. Here, the 




(a), u/t = 0.20t (b), u/t = OAOt (c), and u/t= 0.60* (d). 

same result is obtained as seen from Fig. [3] where the 
panels show constant-energy cuts of x" (<?,<*;) at varying u> 
with fixed U = 2.5t and S = 0.05. 

The discussion above, and the associated results pre- 
sented in FigsOJHJ were all based on a ground state consist- 
ing of a homogeneous d-wave superconductor. Within the 
present mean-field formalism this is the relevant regime 
for clean systems with U < U C 2 — 2.6t. In the opposite 
regime, where U > U C 2, the clean phase is an ordered 
smectic stripe phase which disorder may turn into a ne- 
matic as discussed e.g. in Refs. [45H47] . We do not ex- 
pect a small explicit symmetry breaking term to qualita- 
tively alter this scenario. Instead, we focus on the weak-£7 
limit where disorder may induce a stripe magnetic phase 
as shown in Refs. [38lE9l|4Tl|42] , and study the evolution 
of the nematic spin response to disorder. For this purpose 
we need to investigate how impurities respond to a finite 

Figure Ufa) shows the Ci symmetric pattern of local 
magnetic order induced around such an impurity, which 
we will refer to as a "spin nematogen" , while in Figure 
0Jb) we see that a significant anisotropy is absent in the 
charge sector, as may be expected in weak-coupling the- 
ories of this type. Similar to the bulk results in Fig. [TJ 
the impurity-induced magnetization becomes increasingly 
asymmetric as one approaches V C 2 from below. The ques- 
tion arises how to detect these spin nematogens in the 
cuprates. NMR is clearly sensitive to the local distribution 
of spins, and the lineshapc in other disordered, correlated 
materials has been analyzed in terms of postulated struc- 
tures of this type [32]. In addition, one might expect a 
signature in the local density of states around an impurity; 
the well known fourfold-symmetric impurity resonance as- 
sociated with a potential scatterer in a d-w&ve supercon- 
ductors has been detected by STM [48] in near-optimally 
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Fig. 4: (Color online) Magnetization (a), charge density (b), 
and local density of states (c) plotted in real-space near a single 
strong scatterer (Vi mp = lOOt) with hopping asymmetry So — 
0.05 and U/t = 2.5. In (a-c) we have used N = 35, but show 
only the center 11x11 sites near the impurity for clarity. In (c) 
the LDOS is shown at ui = — 0.05f which is the energy of the 
impurity resonance for the present choice of parameters, (d) 
Constant energy cut of the susceptibility \ '{q,uj) with U/t = 
2.5 and uj/t = O.Ot for the case of disorder-induced nematogens 
corresponding to (a), (e) shows q x (q y ) cuts of x"(q> u ) along 
q y = 7r [black line] (q x = n [red line]) for the same case as 
shown in panel (d). (f) same as (e) but in the homogeneous 
case similar to Fig. [3ja). 



doped samples. In the presence of x — y anisotropy, this 
resonance acquires within our theory a similar anisotropy 
as shown in Fig. @Jc) which should be observable near 
strong scatterers in the underdoped regime. 

The spin nematogens around an impurity arise because 
of a "freezing" of gapped incommensurate spin fluctua- 
tions present in the clean system [3311351135] ■ In agree- 
ment with the x — y anisotropy shown in Fig. BJa), 
the low-energy disorder-induced spin excitations also be- 
come highly anisotropic, as shown in Fig. Eld). Because 
of numerical limitations the spin-susceptibility from in- 
homogeneous real-space configurations is obtained from 
systems with N = 24. Surprisingly, the disorder signif- 
icantly enhances the nematic response as seen by com- 
paring Figs. H£e,f); (e) displays two distinct cuts along 
either q x = 7r or q y = n from Fig. HJd), and reveals a 
much larger anisotropy in the spin response for the dis- 




Fig. 5: (Color online) (a) Static magnetization shown in real- 
space for one of the disordered systems used in (b-d). As 
seen, the frozen magnetization roughly consists of anti-phase- 
coupled antiferromagnetic domains, (b) Constant energy cuts 
of the susceptibility y^'{q,u) with U/t = 2.5 and u/t — O.Ot 
averaged over ten different 24 x 24 systems each containing 
4% impurities of strength Vi = 5t. (c) same as (b) but at 
ui/t — OAOt. (d) shows q x (q y ) cuts of ■x/'{q,uj) from (b) along 
q y = 7r [black line] (q x = n [red line]). 

ordered case (e) compared to the clean case (f). Since a 
gapped spectrum near the Fermi level (caused by super- 
conductivity or pseudo-gap physics) is a prerequisite for 
disorder-induced magnetization, we obtain the interesting 
situation where a disordered superconductor in the pres- 
ence of sub-dominant electronic correlations and a small 
x — y symmetry breaking field, work as a catalyst for ob- 
serving a nematic response in neutron scattering. 

We end this section with a brief discussion of the real- 
istic many-impurity situation where the impurity-induced 
magnetization form a glassy pattern shown in Fig. Efa). 
Panels (b) and (c) in Fig. [5J displays respectively the low- 
and high-energy spin susceptibility, '(q,oj), in the pres- 
ence of a few percent added disorder similar to the exper- 
imental study in e.g. Ref. [49]. The particular concen- 
tration of disorder (and their strength) is not important 
for the present discussion. In agreement with the neutron 
measurements, the low-energy spin response from a collec- 
tion of overlapping spin nematogens peaks at the incom- 
mensurate (commensurate) position along q x (q y ). Again 
this distinct asymmetry is only present in the low-energy 
sector as seen from Fig. [FJc). In Fig. [FJd) we show cuts 
through Fig. HJb) along q y = ir [black line] and q x = it [red 
line] , revealing a semi-quantitative fit to the experimental 
low-energy spin fluctuations obtained in Ref. |17j . 

Conclusions. — We have calculated the spin suscep- 
tibility of a d-wave superconductor with Hubbard correla- 
tions in the presence of small explicit symmetry breaking 
of the underlying lattice. Correlations were found to sig- 
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nificantly enhance the low-energy anisotropy of the spin 
response as the stripe instability is approached. Compar- 
isons with experiment [17l|49] , as well as with work based 
on the same model and recent analysis of the doping de- 
pendence of such effects at strong coupling [50] suggest 
that this effect is responsible for the strongly enhanced 
nematic tendency in the spin response observed in the 
YBCO system as it is underdoped. In addition, we have 
shown here that disorder significantly enhances the ne- 
maticity, via generation of local spin nematogens which 
exhibit an anisotropic low-energy spin response that can 
be significantly enhanced compared to the clean case, and 
is crucial to understand the g-space form of the neutron 
response near (n, n) observed in experiments on strongly 
underdoped, untwinned YBCO samples. Our predictions 
can be tested by studying the effect of Zn substitution 
on the anisotropic spin response, and by analysis of NMR 
and STM spectroscopy of impurity bound states in under- 
doped cuprates. 
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